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Abstract. We consider the recursion operator approach to the soliton equations related to 
the generalized Zakharov-Shabat system on the algebra sl(n, C) in pole gauge both in the 
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1 Introduction 

m 

The theory of nonlinear evolution equations (NLEEs) of soliton type (soliton equations or com- 
pletely integrable equations) has developed considerably in recent decades. Interest in it is still 
' big and there is a wide variety of approaches to these equations. However, some properties are 

T-H fundamental to all approaches and one is that these equations admit the so called Lax repre- 

sentation, namely [i,^] = 0. In the last expression L and A are linear operators on dx, dt 
. . depending also on some functions qaix,t), 1 < a < s (called 'potentials') and a spectral para- 

meter A. Since the Lax equation vl] = must be satisfied identically in A, it is equivalent to 
^ a system (in the case when A depends linearly on dt) of the type 

^ {qa)t = Fa{q,qx,- ■ ■), where q = {qa)i<a<s- (1) 

In most of the approaches the linear problem Lip = (auxiliary linear problem) remains fixed 
and the evolution equations (of a certain form) that can be obtained by changing the opera- 
tor A are considered. The hierarchies of equations we obtain by fixing L are called the nonlinear 
evolution equations (NLEEs), or soliton equations, associated with (or related to) L (or with the 
linear system Lip = 0). The hierarchies usually are named for some of the remarkable equations 
contained in them. The schemes according to which one can calculate the solutions to the soliton 
equations may be very different, but the essential fact is that the Lax representation permits 
one to pass from the original evolution defined by the system of equations (1) to the evolution 
of some spectral data related to the problem Lip = 0. Since finding the spectral data evolution 
usually is not a problem, the principal difficulty is to recover the potentials from the spectral 
data. This process is called the inverse scattering method, which is described in detail in the 
monographs [4, 10]. 
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The generalized Zakharov-Shabat system (GZS system) we see below is one of the best known 
auxiliary linear problems. It can be written as follows 



Here q{x) and J belong to a fixed simple Lie algebra q in some finite-dimensional irreducible 
representation. The element J is regular, that is the kernel of adj (adj(X) = [J, X], X £ q) is 
a Cartan subalgebra f) C g. The potential q{x) belongs to the orthogonal complement t}^ of f) 
with respect to the Killing form 



and therefore q{x) = QaEa where Eq, are the root vectors; A is the root system of g. The 

scalar functions qa{x) (the 'potentials') are defined on M, are complex valued, smooth and tend 
to zero as X — )• ±00. We can assume that they are Schwartz-type functions. The classical 
Zakharov-Shabat system is obtained for g = sl(2, C), J = diag(l, —1). 

Remark 1. We assume that the basic properties of the semisimple Lie algebras (real and 
complex) are known and we do not give definitions of all the concepts related to them. All our 
definitions and normalizations coincide with those in [15]. 

Remark 2. When generalized Zakharov-Shabat systems on different algebras are involved we 
say that we have a generalized Zakharov-Shabat g-system (or generalized Zakharov-Shabat 
on g) to underline the fact that it is on the algebra g. When we work on a fixed algebra its 
symbol is usually omitted. 

Here we may mention also that, in case when the element J is complex, the problem (2) is 
referred as a Caudrey-Beals-Coifman system [2] and only in the case when J is real it is called 
a generalized Zakharov-Shabat system. The reason for the name change is that the spectral 
theory of L is of primary importance for the development of the inverse scattering techniques 
for L and the cases when J is real or complex are quite different from the spectral viewpoint, 
see for example [5, 11] in which the completeness of the so-called adjoint solutions of L when L 
is considered in an arbitrary faithful representation of the algebra g is proved. Referring for 
the details to the above work we simply remind the reader that the adjoint solutions of L are 
functions of the type w = mXm~^ where X is a constant element from g and m is a fundamental 
solution of Lm = 0. Let us denote the orthogonal projector (with respect to the Killing form (3)) 
on f)-*- by ttq. Then, of course, the orthogonal projector on t) will be equal to id — ttq. Further, let 
us put = ttqw and w*^ = (id —tto)w. One of the most important facts from the theory of GZS 
systems is that if a suitable set of adjoint solutions {wi{x,X)) is taken, then, roughly speaking, 
for A on the spectrum of L the functions {wf{x, A)) form a complete set in the space of potentials. 
If one expands the potential over the subset of the adjoint solutions the coefficients are one of the 
possible minimal scattering data sets for L. Thus, passing from the potentials to the scattering 
data can be considered as a sort of Fourier transform, called a generalized Fourier transform. 
For it wf{x, A) play the role the exponents play in the usual Fourier transform. They are called 
generalized exponents or by abuse of language also adjoint solutions. Those familiar with the 
theory in sl(2; C) know that originally they were called 'squares' of the solutions of Lip = 0. 
This interpretation of the inverse scattering transform was given for the first time in [1] and 
after that has been developed in a number of works, see, for example, the monographs [10, 17] 
for complete study of sl(2,C)-case and for comprehensive bibliographies, and [2, 11] for more 
general situations. 

However, since in this article we shall not deal with the spectral properties of L, we shall call 
it a generalized Zakharov-Shabat system in all cases. 



= {id^ + q{x) -XJ)ib = 0. 




{X,Y) = tr(adxady) 



(3) 
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The recursion operators (generating operators, A-operators) are the operators for which the 
functions wf{x, A) are eigenfunctions and therefore for the generahzed Fourier transform they 
play the same role as the differentiation operator in the usual Fourier transform method. For that 
reason recursion operators play central role in the theory of soliton equations - it is a theoretical 
tool which, apart from explicit solutions, can give most of the information about the NLEEs 
[10, 36]. The theory of these operators is an interesting and developing area. Through them one 
may obtain: 

i) the hierarchies of the nonlinear evolution equations solvable through L; 

ii) the conservation laws for these NLEEs; 

iii) the hierarchies of Hamiltonian structures for these NLEEs. 

It is not hard to find that the recursion operators related to L have the form, see [11] or the 
book [10] and the numerous references therein, 

A±iX{x)) = adf (^id,X + Tro[q,X]+iadg£ {id-7roMy),X{y)]dy^ . (4) 

Here, of course, adq{X) = [q, X] and X is a smooth, rapidly decaying function with values in f)-*". 

The name 'recursion operators' has the following origin. Suppose we are looking for the 
NLEEs that have Lax representation [-Z>, ^] =0 with L given in (2) and A of the form 

n 

A = idt + Yl ^" ^ ^' ^" = ^"-1 ^ 

k=0 

Then from the condition [-^^, j4] = we first obtain An-i = adj^[g, A„] and for < k < n — 1 
the recursion relations 

TTo^fc-i = A±(7ro^fc), (id -7ro)^fc = i(id -TTo) / [q,TToAk]{y)dy, 

J itoo 

where A± are as in (4). This leads to the fact that the NLEEs related to L can be written in 
one of the following equivalent forms: 

a) iadj^qt + Al {ad j^[An,q]) =0, (5) 

b) iadfqt + A'l{adf[An,q])=0. 

Remark 3. Strictly speaking, this is not the most general form of the equations solvable 
through L. Considering the right-hand side of the equations of the type ad qt = Fn{q) as 
vector fields, in order to obtain the general form of the NLEEs associated with L one must take 
an arbitrary finite linear combination F of the vector fields Fn with constant coefficients and 
write adj^ qt = F{q). We refer to (5) as the general form of the equations solvable through L 
for the sake of brevity. 

There is another important trend in the theory of the GZS system and consequently for the 
recursion operators related to it. It turns out that this system is closely related to another one, 
called the GZS system in pole gauge (then the system L we introduced is called GZS system in 
canonical gauge). In order to introduce it, denote the group that corresponds to the algebra g 
by G. Then the system we are talking about is the following 



Li> = id^i> -Si> = 0, S gOj 



(6) 
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(with appropriate conditions on S{x) when x i— t- ±c«), where Oj is the orbit of the adjoint 
action of the group G, passing through the element J € g. A gauge transformation of the type 
ip I— >• ipQ^ip = 4' where ■00 is a fundamental solution to the GZS system corresponding to A = 
takes the system L into the system L if we denote S = ipQ Jipo- One can choose different 
fundamental solutions 'i/'o and one will obtain different limiting values for S when x i— >• ±00 but 
usually for ■00 is taken the Jost soliton that satisfies lim ipQ = 1. The GZS system in pole 

x—^—oo 

gauge is used as the auxiliary linear problem to solve the equations that are classical analogues 
of equations describing waves in magnetic chains. For example, in the case of g = sl(2, C) 
(spin 1/2), one of the NLEEs related to L is the Heisenberg ferromagnet equation 

St = j^[S,Sxx], 5Gisu(2). (7) 

In the case of sl(3, C) the linear problem L is related to a classical analog of equations describing 
the dynamics of spin 1 particle chains [3, 31]. 

The theory of NLEEs related to the GZS auxiliary system in canonical gauge (L) is in 
direct connection with the theory of the NLEEs related with the GZS auxiliary system in pole 
gauge (-L). The NLEEs for both systems are in one-to-one correspondence and are called gauge- 
equivalent equations. This beautiful construction was used for the first time in the famous work 
of Zakharov and Takhtadjan [37], in which the gauge-equivalence of two famous equations - the 
Heisenberg ferromagnet equation (7) and the nonlinear Schrodinger equation was proved. 

In fact the constructions for the system L and its gauge equivalent L are in complete analogy. 
Instead of the fixed Cartan subalgebra f) = keradj, we have a 'moving' Cartan subalgebra t}^ = 
kerad5(2.); a 'moving' orthogonal (with respect to the Killing form) complementary space t)^ 
to f)5 etc. We have the corresponding adjoint solutions rh = ipX^~'^, where ip is a fundamental 
solution of Lijj = and X is a constant element in g. If we denote by fh'^ and rh'^ the projections 
of rh on i)^ and f)^ respectively, then the corresponding recursion operators are constructed 
using the fact that the functions must be eigenfunctions for them. The evolution equations 
associated with the system (6) and gauge-equivalent to the equations (5) have the form 

-i ads' St + {A±rA±7rsAn = 0, (8) 

where A.± are the recursion operators for L and vrg is the orthogonal projector on f)^. One can 
see that 

A± = Ad (■00 ^) o A± o Ad(0o), T^s = Ad {ipQ^) onoo Ad{7po), 

where Ad is denotes the adjoint action of the simply connected Lie group G having g as algebra. 

Remark 4. In order to understand why in the hierarchy (8) appears TTsAn one must mention 
that one can prove that for any constant H £ i) 

Ad {^Q^) {adf [q,H]) =7rsH. 

So for the GZS system in pole gauge everything could be reformulated and the only difficulty 
is to calculate all the quantities that are expressed through q and its derivatives through S and 
its derivatives. Though in each particular case the details may be different there is a clear 
procedure to achieving that goal. The procedure was developed in the PhD thesis [30], outlined 
in [12, 14] for the sl(2,C) case, and for more general cases in [13]. In the case sl(3,C) the 
procedure has been carried out in detail in [31]. The theory of the recursion operator for the 
GZS system in canonical gauge in the presence of the Mikhailov type reductions has been also 
considered, see for example [6] for a treatment from the spectral theory viewpoint, or [33] for 
geometric treatment. 
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The theory for the system in pole gauge in presence of Mikhailov type reductions has been 
also the subject of recent research. In [7, 8] the case sl(3, C) in the presence of Z2 x Z2 reduction 
was considered (Gerdjikov-Mikhailov-Valchev system or GMV system). In [34] it was shown 
that the operators found in [8] using classical technique and the technique developed in [16], 
are restrictions of the general position recursion operators for sl(3, C), on certain subspaces of 
functions. In [35] the geometric theory behind the recursion operators for the GMV system was 
presented. In the present article we shall generalize the theory developed in [33, 34, 35] for the 
algebra sl(3,C) to the algebra sl(re, C). Another new feature is that we discuss here in more 
depth the fundamental fields of the Poisson-Nijenhuis structure (P-N structure) related to the 
GMV system and its generalizations. 



2 Recursion operators for the GZS sl(n, C)-system in pole gauge 
2.1 Preliminary results 

Let us consider the GZS system in canonical gauge in the case g = sl(n, C) ~ the algebra of 
all traceless n x n complex matrices. The Cartan subalgebra f) consists of all traceless diagonal 
matrices and the space t)^ of all off-diagonal matrices. The Killing form can be expressed 
through the trace form, we have {X,Y) = 2ntrXY for any X,Y £ sl(n, C), see [15]. Since the 

n 

element J belongs to t) it has the form J = diag(Ai, A2, . . . , A„), where ^ = 0. Next, J is 

k=l 

regular, which means that for any i ^ j we have Aj 7^ Xj. In that case J 'generates' the Cartan 
subalgebra f) in the following sense. Consider the matrices 

J = Jl, J2 = J 1, Jn-l = J" 1. 

n n 

Since tr = for k = 1, 2, . . . , n these matrices belong to \). One can easily show that these 
matrices are linearly independent and hence generate the Cartan subalgebra. The same can also 
be deduced if one calculates the determinant of the Gram matrix 

/ (Ji,Ji) (Ji,J2) ... (Jl,J„-l) \ 

(J2,Jl) {J2,J2) ■■■ {J2,Jn-l) j-g-j 

XiJu-liJl) {Jn-l,J2) ■■■ {Jn-1, Jn-l)/ 

Using that {X, Y) = 2n tr XY we obtain 

detG = 2"-in"-2 J|(Ai - Xjf = T"-^ nJ'-'^ D'^ . 

i<j 

The determinant of the Gram matrix (9) enters into the explicit calculation of the recursion 
operators (because we need the inverse of G) and as we see it is proportional to the square of 
the Vandermonde determinant D constructed from the values Ai, A2, . . . , A^. This observation 
has been made for the case of the algebra sl(3, C), see (18). We see now that it is a general 
result for sl(n, C). 

Since the Killing form for a given simple Lie algebra is invariant under the adjoint action of 
the corresponding group, we get that the Gram matrix {{Si, Sj))i<ij<:n-i where 

5 = 5i =Ad(V^^)Jl, S2 = Ad{%b^^)j2, Sn-l=Ad{lP^^)Jn-l 

coincides with G. Consequently its determinant is different from zero. 
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Another observation that we want to make is that J and S satisfy the equation 
n 

Yl{X-\^l)=0. 



i=l 



(Since Aj are the eigenfunctions of J and S the above is just the Cayley-Hamilton theorem.) 
This equation can be written in the form 

X" - aiX"~i + a2X"-2 + . . . + (-l)"a„l = 0, 

where the coefficients are homogeneous polynomials of degree s in Ai, A2, • • • , A,i. If necessary, 
one can express them in terms the symmetric polynomials {Js,Js), i.e., polynomials in Aj, 
1 < i < n, using the Newton formulae but it is not needed for our purposes. We only note that 
ai = tr J = 0. 

At the end of our preparations let us denote Sg-^x = dxSg- Then we have: 
Proposition 1. The matrices Si-x, S2-X-, ■ ■ ■ , Sn-i;x belong to 

Proof. We need to prove that for all 1 < k,l < n — 1 we have {Sk,Si-x) = 0. Since Sk = 
— (tr J*^/n)l we must prove that tT{S''{S^)x) = 0. Using the properties of the trace we see 
that this is equivalent to tic{S^~^''~^Sx) = 0. On the other hand, for any integer m > we have 
tr5"* = tr J™ = const. Therefore, dxtr S"^ = mtr S"^"^ Sx = and our result follows. ■ 

2.2 Calculation of the recursion operators 

Now we pass to the calculation of the recursion operator(s) A± for the GZS system in pole gauge 
on sl(n, C). We shall use the equation 

- X[S,w] = 0, 

which is satisfied by every function of the type w = tijA{'ijj)~^ , where A is a constant matrix 
and ip is a fundamental solution of (6). We have 

(Note that kerad^ = i)s, the space ijg is its orthogonal space with respect to the Killing form 
and by upper indices 'h' and 'a' we denote projections onto these spaces.) We have seen that 
the matrices {Si}^~^ (where Si = S) span f)^. Therefore 



n-l 

= ^afc(x)5fc, 

i=l 

where afc(x) are coefficient functions. We have 



n-l 

W" 

k=l 



A[S,?i^] = 0, 



or in other words, 

71—1 71—1 



idxw'' + i Y{ak)xSk + i ^ akSk;x - A[5, u)^] = 0. (10) 



k=l k=l 
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In order to find the coefficients we calculate the inner product of the left hand side of (10) 
with Sj. Then, taking into account Proposition 1, we arrive at the following system 



n-l 



k=l 

where Gik = {Si, S^) = {Ji, Jk) are the entries of the matrix G introduced earlier. Equivalently, 

n-l 

akx = {ak)x = ^ {G'^)^^{w'', Ss-x)- 

k=l 

Assuming that for the eigenfunctions of we have lim = and for the eigenfunctions 
of A- we have lim = we find that 

x— oo 

n-l 



ak = Y.iG-'),d-Hw\Ss 



k=l 

where stands for one of the two operators 



■dy, / 

,-oo J —( 



■dy. 

Consequently, inserting the functions Uk into (10), we obtain A±{'w"') = Xw"" where 

A±(Z) =iad^i7r5 \^xZ+ |] {G~^) ^d'^Z , Ss;x)Sk;}l , (11) 
or equivalently, 



k,s=l 



^±iZ) = iadc TTs < 



n-l 



d^Z- ^ {G-')^d~'{Z,,,Ss)Sk;. \ . (12) 



fc,s=l 



The above operators act on functions Z{x) that are smooth, rapidly decaying and such that 
Z{x) G ^s{x)- '^^^ formulae (11), (12) give us the recursion operators but they can be written 
in more concise form if we introduce: 

• The row vectors 

• The column vectors 



{Z, Sx) — [{Z,Sl-x), {Z, S2;x), ■ ■ ■ ,{Z, Sn-l;x)) , 
{Zx, S) = {{Zx, Si), {Zx, S2), . . . , {Zx, Sn-l)) ■ 



Then the recursion operators acquire the form 

A±{Z) =iads'7Ts{dxZ + SxG-'d-'{Z,Sx)} (13) 
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or 

~k±{Z)=ia.ds^T:s{d,.Z -S^G-^d-\Z^,S)]. (14) 

It is easy to check now that the recursion operators for the GZS systems on the algebras 
sl(2,C), sl(3,C) (see [12, 31, 32]) are obtained from the general expressions (13), (14) as partic- 
ular cases. 

The last thing that remains to be done is to express the operator ad^^ through S. For this 
note that, if all the eigenvalues of J are different, the operators adj and adj^ considered on f)-*- 
(ad^ and ad^^ considered on f)^ respectively) are simple and have common eigenvectors. Then 
we can apply the following proposition which is actually the spectral decomposition theorem for 
a given simple matrix A, see [18]. 

Proposition 2. Let A he a simple matrix with eigenvalues Ai,A2,...,Am over M or C. Let 
/xi, /i2, . . . , ^ra be arbitrary numbers. Let us define the polynomial 



f{x) = Y,i^MX), 



k=l 

where Ik are the Lagrange interpolation polynomials 
(A - Xi) 



h{x) = n 



Then the matrix f{A) has as eigenvalues /ii,/X2, . . . and the same eigenvectors as A and the 
polynomial /(A) is the polynomial of minimal degree having that property. 

Remark 5. It is not difficult to see that lk{A) is the projector onto the subspace corresponding 
to the eigenvalue A^ in the splitting of the space into eigenspaces of the matrix A. 

Remark 6. In case the matrix A is not simple one again can produce a polynomial /(A) of 
minimal degree having the property stated in the Proposition 2 though its construction is more 
complicated, see [18]. In case just adj^ is needed, one can also use the following procedure: 

1. The minimal polynomial m(A) of adj (on the whole algebra) is a product m(A) = Ami (A) 
and A and mi (A) are co-prime. 

2. The algebra splits into direct sum of invariant subspaces ker adj and imadj (because adj 
is skew-symmetric) and the minimal polynomials of adj on these spaces are A and mi (A) 
respectively. 

3. On imadj the operator adj is invertible and one can find adj^ as polynomial in adj 
multiplying the equation mi (ad j) = by adj^. 

Note that the polynomial g{X) such that adj^ = ^(adj) will have now degree deg(mi) — 1. 

In the case of an arbitrary semisimple Lie algebra, assuming that all the values a (J), a G A 
are different, the operators adj and adj^ have eigenvalues a{J) and l/a(J), a G A respectively 
and eigenvectors Eq,, a G A. Then adj^ is equal to /(adj), where /(A) is the polynomial 



In our case the roots are aij = Si — Sj, i ^ j where on diagonal matrices A with diagonal 
elements we have ei{A) = ai. So aij{J) = Aj — Xj and the expression for ad^^ can be 
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written easily, we shall do it a little further. What we obtained already allows us to make 
an important observation. As the polynomial Z(A) from (15) is of the form A/o(A), where /q is 
another polynomial, then we have /(adj) = /o(a'dj)adj. Therefore we have /(adj)7ro = l{adj). 
In the same way 

l{ads) = &d^^, l{ads)TTs = l{ads), tts = adj^ ads • (16) 

Therefore, if one assumes that ad^^ is given by a polynomial in ads, writing the projector vrs 
in the expression for the recursion operators (13), (14) is redundant. 

Continuing our discussion about ad^^, as we already remarked, in the case when we need 
only to express ad^^ through ads, instead of Proposition 2 we can use the minimal polynomial 
for ads (restricted to f}-*- of course). Indeed, in the case of an arbitrary semisimple Lie algebra g 
and regular J let us assume that all the values q( J) are different. Then the minimal polynomial 
for ad J on f)-*- has the form 

m(A)= J] (A2-a(J)2) =^afcA'', = 1. 

aeA+ k=0 

Here P = \ dim(f)-'-) = ^(dimg — dimf)) and, because J is regular, 

«o = (-1)^ n "('^)' ^ 0- 

Since rn-(ads) = we obtain 

ad^"*^ = adsi?(ads), t^s = ad^^ ads = ads-R(ads) ads, (17) 
where R{\) is the polynomial 

R{\) = -— (aiA^ + 02 A^ + • • • + apA^P) . 

Both expressions (16) and (17) for ad^^ give the same result in the case where all a(J) are 
different. Indeed, both polynomials A^(A) and Ai?(A) are monic and when A = a(J), a G A 
give l/a(J). Since they are of degree 2p — 1 these polynomials coincide. In particular, in the 
case of the algebra sl(n, C) we get 




Unfortunately, the expressions for ad^^ become very complicated for big n, hampering the 
possible applications. A simplification can be obtained for some particular choices of J, for 
which the minimal polynomial of adj on i)^ has smaller degree than in the case of general J. 
Here are the expressions for ad^^ used up to now in the literature: 

1. Classical Zakharov-Shabat system, g = sl(2,C), J = diag(l, —1) [12, 14] 
1 
4 



ad_5 = - ads • 



GZS system on g = sl(3,C) in general position, J = diag(Ai, A2, A3), Aj 7^ 0. The last 
condition ensures that all a(J)'s are different [31, 32] 

ad.i = /(ads), 1{X) = ^(\^ -Ic^ , (18) 



where 

C2 = Xl + Xl + Ai, D = (Ai - A2)(A2 - A3)(Ai - A3). 
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3. GZS type system on g = sl(3, C) with Z2 x Z2 reduction (GMV system), J = diag(— 1, 0, 1) 
[8, 32] 

adsi = /(ad5), /(A) = -J(A3-5A). 

In this case the fact that some of the eigenvalues of adj on f)-*- (and consequently of ads 
on f);^) are not simple leads to a decrease in the degree of the polynomial /(A). 

3 Geometric interpretation 

Fixing the element J for the GZS g-system in pole gauge, the smooth function <S'(a;) with 
domain M, see (6), is not subject to any restrictions except that S{x) G Oj and S{x) tends fast 
enough to some constant values when x 1— )• ±00. First let us consider the even more general case 
when S{x) is smooth, takes values in q and when x — )• ±00 tends fast enough to constant values. 
The functions of this type form an infinite-dimensional manifold which we shall denote by Ai. 
Then it is reasonable to assume that the tangent space Ts{M) at S consists of all the smooth 
functions X : M 1— )• g that tend to zero fast enough when x 1— >• ±00. We denote that space 
by 5^(0). We shall also assume that the 'dual space' Tg{M) is equal to 5^(0) and if a G T^(A^), 
X G Ts{M) then 

/+00 
{a{x),X{x))dx, 
-00 

where ( , ) is the Killing form of g. 

Remark 7. In other words, we identify Tg{Ai) and Ts{Ai) using the bi-linear form (( , )). We 
do not want to make the definitions more precise, since we speak rather about geometric picture 
than about precise results. Such results could be obtained after a profound study of the spectral 
theory of L and L. In particular, we have put the dual space in quotation marks because it is 
clearly not equal to the space of continuous linear functionals on ^{q). We also emphasize that 
when we speak about 'allowed' functionals on we mean that ^ G Tg{M) ~ did)- 

Now we want to introduce some facts. The first fact is that since we identify Tg{A4) and 
Ts{Ai) the operators 

a^P{X) = [d^a, Q{a) = adsia), S e M, a G r|(A^) (19) 

can be interpreted as Poisson tensors on the manifold M. This is well known, see for exam- 
ple [10], where the issue has been discussed in detail and the relevant references are given. One 
can also verify directly that if Hi, H2 are two functions (allowed functionals) on the manifold 
of potentials Ai then 

{HuH2}p = {{^^dj-^)) ' {Hi,H2}q = [S, ^j^^ 

are Poisson brackets. It is also known from the general theory that these Poisson tensors are 
compatible [10, Chapter 15]. In other words P + Q is also a Poisson tensor. Note that the 
tensor Q is the canonical Kirillov tensor which acquires this form because the algebra is simple 
and consequently the coadjoint and adjoint representations are equivalent. 

Now let Oj be the orbit of the adjoint action of G passing through J. Let us consider 
the set of smooth functions / : M 1— )• Oj such that when x — )• ±00 they tend fast enough to 
constant values. The set of these functions is denoted by M and clearly can be considered as 
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a submanifold of 7W. U S £ M the tangent space Ts{M) consists of all smooth functions X 
that vanish fast enough when x i— )• ±00 and are such that X(x) G Ts(^x){(^ j)- (Recall that Oj 
is a smooth manifold in the classical sense). We again assume that T^{J\f) ~ Tsi^M) and that 
these spaces are identified via (( , )). 

We can try now to restrict the Poisson tensors P and Q from the manifold to the mani- 
fold M. The problem how to restrict a Poisson tensor on a submanifold has been solved in 
principle [24], see also [27, 28]. We shall use a simplified version of these results (see [22, 23]) 
and we shall call it the first restriction theorem: 

Theorem 1. Let M he a Poisson manifold with Poisson tensor P and M <Z Ai he a submani- 
fold. Let us denote hy ] the inclusion map of M into M, hy Xp{M)m the suhspace of covectors 
a G T^{M) such that 

Pm{a) e dim{Tm{M)) = im(djm), meM 

{where im denotes the image), and by T-^{A4)m the set of all covectors at m ^ M. vanishing 
on the suhspace im(djm); m G {also called the annihilator 0/ im(djm) in T^(A^)). Let the 
following relations hold: 

X*p{M)m + T^{M)m = T*^{M), meM, 
X*p{M)mnT^{M)mCkej:{P^), m £ M. 

Then there exists unique Poisson tensor P on M, 'yrelated with P, that is 

Pm = djm o Pm o {dim)* , m £ M. 

The proof of the theorem is constructive. One takes /3 G T^{M.), then represents (j*/3)m 
as ai + 02 where ai G Xp{A4)m, 0L2 G T^{M)m and puts Pm{P) = -Fm(cti) (we identify m 
and i{m) here). 

The restriction we present below has been carried out in various works in the simplest case 
= sl(2, C), see for example [23]. We do it now in the case q = sl(n, C). Restricting the Poisson 
tensor Q is easy, one readily gets that the restriction Q is given by the same formula as before 

Q{a) = ads{a), SeM, a G r|(AA). (20) 

The tensor P is a little harder to restrict. Let us introduce some notation and facts first. Since J 
is a regular element from the Cartan subalgebra t), each element S from the orbit Oj is regular. 
Therefore f)5 = kerad5 is a Cartan subalgebra of sl(n, C) and we have 

sl{n,C) = \js{x)ei)U^)- 

If X G Ts{M) then X{x) G i)s{x) and X vanishes rapidly when x 1— )• ±00. We shall denote the 
set of these functions by 5^(f)5) so X G ^^((l^) (this means a little more than simply X G f)^). 
Using the same logic, for X G i?(f)5) we write ads{X) G d{hs) which means that the function 
ads(x) ^{x) belongs to T?(f)5). Now we are in a position to perform the restriction of P on Af. 
For S £ M we have 

X*p{Af)s = {a : idxa G ^{[)^)}, T^{Ar)s = {« : {{a,X)) = 0, X € d{i)i)}. 

We see that T-^{Af)s is the set of smooth functions a such that a G {)s and such that they 
vanish fast enough when x i-> ±00. We shall denote this space by 5(^)5)- We introduce also 
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the space of functions X E 1)3 which tend rapidly to some constant values when x ±oo and 
denote this space by d{i)s)o- Clearly, since Sk, k = 1,2, . . . ,n — 1 span i)s, we have 



n-1 



d{hs)o = < : = '^^ak{x)Sk{x), ak{x) smooth, ak{x) tend to const as a; i->- ±00 > , 

I k=l J 

^{^)s) = Ix : X = ^ak{x)Sk{x), ak{x) smooth, lim afc(x) = o\ . 



k=l 



In what follows we shall adopt matrix notation and shall denote by A(x) the column with 
components ak{x). Then ^p(AA)5 H T-^{J\f)s consists of the elements 

n-1 

a = XI (ik{x)Sk{x) = SA, 
k=l 

with ak{x) vanishing at infinity and such that idxOt G ^{^g)- But 
\dxOi = iSa;A + iS^^, 

so we must have = Q and as{x) are constants. But as{x) must also vanish at infinity so we 
see that as = 0. Consequently Xp{N)s n T^{N)s = {0} C ker Pg. 

Let us take now arbitrary a G T*{J\f)s- We want to represent it as ai + q;2, ai G X*{N)s, 
a2 G T-^{M)s- First of all, a2 = SB{x) where B{x) a column with components bs{x) - scalar 
functions that vanish at infinity. This means that 

id^a = id^rOii + iS^B + iSB^, 
where idxai G d{^s)- But then we have 

{dxa,S{x)) = GBx, 
where G is the Gram matrix we introduced earlier. Therefore 

B = G-^d-\dxa{x),S{x)). 

Remark 8. In the theory of recursion operators when one calculates the hierarchies of NLEEs 
or the conservation laws the expressions on which the operator acts are total derivatives. 
Thus the same results will be obtained choosing for d^^ any of the following operators 



• dy, / • dy, - / • dy + / • dy 

-00 J +00 ^ \J—oo J +00 J 



However, more frequently one uses the third expression when one writes the corresponding 
Poisson tensors in order to make them explicitly skew-symmetric. 

Returning to our task, let us put 

a = ai-|-a2, a\ = a — a2, a2 = ^G~^ {pxOL{x) ,^{x)) , 

where a\ and lie in the spaces X*{M)s and T^{M)s = 5(f)s) respectively. We note also that 
djj^ = 7r5(/3). Thus the conditions of the first restriction theorem are fulfilled and if /5 G Tg{J\f) 
the restriction P of P on J\f has the form 

Pi/3) = msdxP - iSxd-\dxPix), Six)). 
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The tensor Q is invertible on AA, so one can construct a Nijenhuis tensor N = P o ad^^ 

N{X) = i7rsd^{ads^X)-iS,G-^d-^{d^{a.ds^X),S{x)), X G ^(f)^)- (21) 

Taking into account that {adg^{X), S) = 0, the above can be cast into the equivalent form 

N{X) = iTTsd.iads^ X) + iS.,G-^d~\ads' X, S,), Xg d{i)s)- (22) 

From the general theory of compatible Poisson tensors it now follows that 

Theorem 2. The Poisson tensor field Q (20) and the Nijenhuis tensor field N given by (21), (22) 
endow the manifold M with a P-N structure. 

The final step is to calculate the dual of the tensor A'^ with respect to the pairing (( , )). 
A quick calculation, taking into account that ad^ is skew-symmetric with respect to the Killing 
form, gives 

N*{a) =ieids' [Trsd,a + S,G-\a,S,)] , a G i?(f)^), 
or equivalently, 

N*{a) = iadgi [Ksd^a - S.,G-^d-\a^,S)] , a e d{t)s)- 

But these are the recursion operators A-t for the GZS system in pole gauge from (13), (14) and 
our results confirm the general fact that the recursion operators and the Nijenhuis tensors are 
dual objects. 

Now, according to the general theory of recursion operators, see [10], the NLEEs related to 
the system L have the form 

-iad^i — + (A±)"A±^5^ = 0, (23) 

where H is an element of the Cartan subalgebra of sl(n,C). As discussed earlier, its can be 
expressed as a polynomial in ad^. Then (23) gives the hierarchies of NLEEs related to the 
GZS system in pole gauge as hierarchies of equations gauge-equivalent to hierarchies related 
to the GZS system in canonical gauge. If, however, one is not interested in finding the pairs 
of gauge-equivalent equations, but wants just to find the NLEEs, one can proceed as follows. 
Recall that recursion operators also produce the hierarchy of Lax pairs. In fact, if a NLEE has 
Lax representation [L, ^] = with L in general position and A has the form 



k=0 

then idxAk £ Aq = const. Using a gauge transformation depending only on t one can ensure 
that Aq = and then the coefficients Ak for k = 1, 2, . . . , n — 1 may be calculated with the help 
of the recursion operator in the following way 

Al_, = A^Al, At = iSG-'d^\Al,S,), = (id-7r5)i„ if = ttsA,. 

For the coefficient An-i one has that idxAn = [S,An-i] and, since An € i}s, there are n — 1 
scalar functions ak{x) forming a column vector a{x) such that An = S{x)a{x). This gives 

iSxa{x) + iSaxix) G f)^, 

and therefore = so ak{x) are constants. Thus An-i = iad^^(Sa) and the hierarchy of the 
NLEEs related to L is 



dS 

ad^^ — + (A±)"(ad^^(S:;;a)) =0, a = (01,02, . . . ,a„_i)*. 
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4 Reductions 

4.1 Algebraic aspects 

As a matter of fact the situation most interesting for applications is when we do not have 
the GZS system in pole gauge in the general case, but when additional restrictions are made. 
For example, in the Heisenberg ferromagnet equation we require that S = (where f stands 
for Hermitian conjugation). Similarly, for systems describing spin chain dynamics that we 
mentioned in the case of sl(3, C), we also require = S. All this means that we require for 
n = 2,3 that S belongs to isu(n). The algebra su(n) is a real form of sl(n,C) with respect to 
complex conjugation a{X) = —X^ and 

sl(n, C) = su(n) © isu(n), [su(n), su(n)] C su(n), 
[su(n), isu(n)] C isu(n), [isu(?7,), isu(n)] C su(n). 

The space i su(n) is the space of n x n Hermitian matrices. Of course, it is not a real Lie algebra 
since, li X,Y G isu(n), then [^,1"] G su(n) and hence i[X, y] G isu(n). 
Next, since the Cartan subalgebra [} is invariant under a, it also splits 

f) = (f) n su(n)) © ([) n isu(n)). 

The first space consists of diagonal matrices with purely imaginary entries, while the second 
consists of diagonal matrices with real entries. If we want S to be Hermitian, then it is natural 
to assume that J is real. In fact the reduction one obtains in this way is most effectively treated 
by the notion of Mikhailov's reduction group, see [19, 25, 26]. According to that concept, in 
order to perform a reduction we must have a group Go acting on the fundamental solutions of 
the linear problems L and A in the Lax representation. Assume we have the GZS system in pole 
gauge and we take the group generated by one element gi, acting on the fundamental solutions 
as 

5iWx;A)) = [^"i(x;A*)]t, (24) 

where * stands for complex conjugation. Since gf = id, the group generated by gi is isomorphic 
to Z2 and we have Z2-reduction. The invariance of the set of fundamental solutions under the 
action (24) means that o"(5) = —S (that is S" S isu(n)) and that S belongs to the orbit of SU(n) 
passing through J G isu(n). Let us denote this orbit by Oj{S\J{n)). Thus J must be real. 
Moreover, if the operator A has the form 

N 

A = idt + Y,^kX\ 

k=0 

then, if the set of fundamental common solutions ip of (6) and Atp = is Go-invariant, one must 
have also cr(^fc) = A/^ where a is the complex conjugation introduced earlier, that is G i su(n). 
One may note that the generating operators are the same as before simply we may assume that 
all matrices are in isu(n). 

One can have also another complex conjugation of sl(n,C). Denote it by r. Then one can 
construct another Z2 reduction. One can, for example, take t{X) = X* and g acting as 

g{ijj{x;X)) = {tP{x;X*r. 

Then one obtains that S* = —S, A'^ = —A'^, so introducing = iS and A^ = iA^ one sees that 
(after canceling i) the L, A pair is real, with matrices belonging to sl(n,M). 
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It is possible to have restrictions defined by two complex conjugations. This can be achieved 
in the following way. Suppose r is another complex conjugation, commuting with a. Then as it 
is well known h = ar = ra is an involutive automorphism of the algebra sl(re, C). In fact doing 
the things the other way round is easier - to find an involutive automorphism h commuting 
with a and then put t = ha = ah. One can take, for example, h{X) = JkXJk, where Jk is 
a diagonal matrix 

diag(l,l,...,l,-l,-l,...,-l). 



K times n—K times 

and one can put t = ha = ah. Then one can consider the group Go with generating elements gi 
(as in (24)) and where 

52(V^(x; A)) = h{i){x- -A)) = Jk{^{x- -\))Jk- (25) 

Since gig2 = g29i, 91=92 — id the reduction group Gq is isomorphic to Z2 x so the 
invariance of the fundamental solutions with respect to the actions (24) and (25) defines Z2 x Z2 
reduction. In [8, 9] the above reduction has been applied in the case of the algebra sl(3, C) 
with Ji = diag(l, —1,-1) (here K = \). A similar reduction group but this time related to the 
algebra so(5) has been considered recently in [29]. 

Returning to the general case, note that instead of the generators gi, 52 we could use 171, gig2 
and 

gig2{4>{x-\)) = JK[iJ-\x;-X*))]^JK. 
The invariance with respect to the group generated by gi, 52 means that 

a{S) = -S, h{S) = -S, a{As) = -As, h{A2k) = ^2fc, h{A2k-i) = -^2fc-i, 
or equivalently 

a{S) = -S, t{S) = S, a{As) = -As, r(^2fe) = -^2fc, ^(^2^.1) = ^2fc-i- 

The complex conjugation r also splits the algebra sl(n, C) 

sl(n, C) = s © L5, [5,5] Cs, [5, is] C is, [is,is]cs. 

(On s the map r is equal to id, on is to — id.) We prefer to work with the involutive automor- 
phism h instead of r. Then the algebra splits into two invariant subspaces for h and these spaces 
are orthogonal with respect to the Killing form. We have 

U = {X :hiX) = {-irX}, e = 0, 1, sl(n, C) = fo © fi- 

A calculation shows that the subalgebra fo consists of matrices Uq having block form 

Uo = , where tr yl + tr S = 0, 

and the diagonal blocks A and B have dimensions K x K and (n — K) x (n — K) respectively. 
In terms of the same type of block matrices the space fi consists of matrices of the type 



C 
D 
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Since a and h commute the spaces of the real and purely imaginary elements for u - namely 
su(n) and isu(n), are split by h into two subspaces 

su(n) = (fo n su(n)) (fi n su(n)), isu(n) = (fo) © (fi). 

As we shall see S{x) G fi so the space fi is of particular interest to us. As is easily seen, this 
space consists of matrices of the form 

^^ = (ct o)' 

while fo n i su(n) consists of matrices of the form 

^°"(o s)' ^'^=A B^=B, tr^ + tr5 = 0. 

Therefore the real vector space fi fl i su(n) can be considered as isomorphic to the quotient space 
su(n)/(s(u(K) X v.{n-K))). 

Below, since all the matrices S'fc, S]^-^x belong to isu(n) (that is, they are Hermitian), we shall 
not write explicitly that they belong to isu(n) and shall concern ourselves only with whether 
they belong to fo, fi, 1)5 and f)^. This makes the formulas simpler and is possible due to the 
facts that restricting to isu(n) does not change the form of the recursion operators and the 
restriction of the Killing form to i su(n) is again nondegeneratc. 

Now, since /i(S') = —S we have ad^o/i = —ho adg and the Cartan subalgebra f)5 = keradg' 
also splits into two subspaces orthogonal to each other 

[3s = %5®'f)s, f3s = f)snfo, 'f)s = D5nfi. 

Since orthogonality with respect to the Killing form is preserved by /i, the space t)^ is also 
invariant under h and 

1)5 =°f)^©'f)5, %5 = ()5nfo, ^f)^ = f3^nfi. 

Consider now S = Si, S2, ■ ■ ■ , Sn-i, the basis of f)5 we introduced earlier. Prom the above it 
follows that 

Sie'i)s, S2e%s, S^e'tjs, 
and one can write 

Sk &%s^ ^ even, Sk &^i)sy ^ odd. 
Consequently, 

Sk;x e k even, Sk;x e ^i)5, k odd. 

In addition, since fo and fi are orthogonal with respect to the Killing form, we have that 

{S2s, S2k-l) = 0. 

Finally, one immediately sees that 
Let us introduce the following spaces: 
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1- ^^(^f)^) consists of all smooth, rapidly decaying functions X{x) on the line such that 
X{x) G ^^s(x) write simply X G ^i)s)- 

2. di^^s) consists of all smooth, rapidly decaying functions X{x) on the line such that 
X{x) G ^^s(x) write simply X G ^hs)- 

Naturally, 

where the spaces are orthogonal with respect to the form (( • )). 

Remark 9. Note that all our matrices are also elements from isu(n) so strictly speaking Sk are 
elements from t}5.nisu(n) or f)^nisu(n) and Sk-x are elements from '^fj^r nisu(n) or -"^fj^r nisu(n) 
but we agreed not to write isu(n) as this will make the notation even more complicated. 

After these preliminaries, assuming that all the quantities are as above we have 

Proposition 3. The recursion operators A± interchange the spaces ^{^i)s) CLn-d dChs) in- the 
sense that 

Proof. In order to make the calculation easier, let us introduce the rows of elements: 

Si = (^l, 5*3, ... , S2k-l), Sl-x = (Sl-x, Ss-x, ■ ■ ■ , S2k-l;x), 

where 2fc — 1 is the largest odd number less then or equal to n and 

So = {S2, 5*4, . . . , 52s), So;a; = {S2-x, S^-x, • • • , S2s;x), 

where 2s is the largest even number less then or equal to n. Let us also introduce the Gram 
matrices 



/ {82,82) {82,84) 
{84,82) {84,84) 



\{82s,82) {823,84) .. 

( {81,81) {81,83) 
{83,81) {83,83) 



{82, 82s) \ 

{84, 82s) 



{82s, 82s)/ 

{8i,82k-i) \ 

{82,82k-l 



(27) 



(28) 



\{82k-l,8l) {82k-l,83) ... {82k-l, 82k-l) / 



Since {Jm,Jn) = {8m, 8n) these matrices have constant entries. With the help of the matri- 
ces (27), (28) the recursion operators can be written into the form 

A±(Z) = iad^i 7:s{dxZ - SQ,x('G)-^d-\Zx, So) - Si,x{^G)-^d-\Zx, Si)} 



or 



A±(Z) = iad^i 7:s{dxZ + S^.x^" G)-^ d'^Z , So,^) + Si-^i' Gy'd'^Z, Si;,)}. 
Assume that Z G dChs)- Then since {Z, So) = {Z, Si) = and {Z, Si-x) = we get that 



A±(Z) = iad^^ TTsidxZ - So-xCGy d-'{Zx, So)} = Al{Z), 
A±(Z) = iad^i 7Ts{dxZ + So-xCG)-'d-\Z, Sq;,)} = A^ (Z). 



(29) 
(30) 
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Recalling (26) we get that h.±{Z) G 5^(^f)5). The above also means that trough the expres- 
sions (29), (30) we define operators - the restrictions of K± on S^C^f}^)- Similarly, if 
Z G ^("^f)^) we have 

A±(Z) =iad^i7r5RZ-Si;,(^G)-ia-i(Z,,Si)} = Ai(Z), (31) 
A±(Z) = iad^i ^s{d.Z + ^i,^{^Gr^d~\Z, Si;.)} = k\{Z), (32) 

and we get that A.±{Z) G 5^('^f)5). In the same way as above (31), (32) define operators Aj_ - 
the restrictions of A-t on 'Si^^^)- other words we have 

A±(5(%^)) e 5(^f)^), A±(5(if)^)) G ^(%^). ■ 

In similar situations (when we have Z2 reductions) the operators A^., K\. are considered 
as factorizing the recursion operator. In some sense this is true, since if one considers (A-t)^ 
acting on S^Cfls) we can write it as A5j_A^. On the other hand (A-t)^ acting on ^^(^f)^) can be 
written as Aj.A^j_. We think that is more accurate to treat the operators A^, A3j_ as restrictions 
of the operators A±. The geometric picture we are going to produce below also supports this 
viewpoint. 



4.2 Geometric aspects 

The geometric situation in the presence of reductions is also interesting. The point is that the 
canonical Poisson structure Qs = ads simply trivializes, apparently destroying the geometric 
interpretation given in the case of the GZS pole gauge system in general position. Indeed, under 
the restrictions considered in this section we first note that the space on which the 'point' S{x) 
takes its values is iC'j(SU(n)) where J is real (and regular). As remarked already this simply 
makes all the matrices Hermitian and everything remains as it was. However, imposing the 
second Z2 reduction (the one defined by the automorphism h) means that S{x) belongs to the 
space of matrices taking values in fi and such that they converge rapidly to some constant values 
when X 1— t- iLcxd. So for the manifold of potentials Q we have: 

Q = |5 : 5(x) G (iC'j(SU(n)))nfi, lim 5 = 5(±oo), fast enough). (33) 

Then the tangent space Ts{Q) to the manifold (33) at the point 5 is the space (for the 

sake of brevity we again 'forget' to mention that the tangent vectors must be also elements of 
isu(n)). But because ad^ interchanges 5(^^)5) and ^i^^g) and they are orthogonal with respect 
to (( , )) for Zi,Z2 G Ts[Q) we have that {{Zi^ [S, Z2])) = 0, so indeed the tensor Q becomes 
trivial. Let us see what happens with the Nijenhuis tensor. Writing everything with the notation 
we introduced in this section we have 

Af(X) = i7r5a.(ad5^X) -iSo;z(°G)-i9-i(a.(ad^i A),So(x)) 

= msd.{B.ds^ X) + iSo;z(°G)-i9-i(ad5i X, Sq;.) = iVo(^) for X G 5(%^), 

N{X) = msd^{B.ds^X)-iSi,^{^Gr^d-\d^{s.d-s^X),^i{x)) 

= msd:,{B.d-s^ X) + iSi;,(iG)-i5-i(ad^i X, Si;,) = Ni{X) for X G ^('f)^)- 

Again we can assume that the right-hand sides of the above equations define the operators A'^i, 
and 
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An immediate calculation shows that for the conjugates of the operators N^, Ni with respect to 
(( , )) we have 

A^o = A±, ^1 = A±- 
Remark 10. One should bear in mind that that enters in the expressions for Nq, Ni is 
• dy or as / • dy. 

-OO J +00 

One can easily see that N{Ts{Q)) does not belong to Ts{Q), so the restriction of A on Q 
cannot be a Nijenhuis tensor on Q. However, N'^idC'^s)) ^ di^hs)^ ^^^^ 

N\TsiQ)) C TsiQ), 

so N'^ becomes a natural candidate. Indeed, let us recall the following facts from the theory of 
P-N manifolds, see [10, 20, 21, 22]: 

Theorem 3. IfAi is a P-N manifold endowed with Poisson structure P and Nijenhuis tensor N, 
then for k = 1,2, . . . each pair {N^P = P{N*)'' , N^) also endows M with a P-N structure. 

Theorem 4. Let Ai be a P-N manifold endowed with Poisson tensor P and Nijenhuis tensor N. 
Let Ai C Ai be a submanifold of M and suppose that we have: 

i) P allows a restriction P on M such that if ] : Ai ^ Ai is the inclusion map then P is 
yrelated with P, that is Pm = djm ° Pm ° (djm)*) rn G Ai. 

a) The tangent spaces of Ai, considered as subspaces of the tangent spaces of Ai are invariant 
under N , so that N allows a natural restriction N to Ai, that is N is ^-related with N. 

Then (P, N) endow Ai with a P-N structure. 

We call the above theorem the second restriction theorem. 

In view of what we have already, we need only find the restriction P' of P on Q and then P' 
and the restriction of A^ will endow Q with a P-N structure. Thus we have the following 
candidates for restriction - the Poisson tensor P = N o ad^ = adg oN* and the Nijenhuis 
tensor N"^ (or N~'^). Let us take P and try to restrict it. 

We want to apply the first restriction theorem. X*{P)s consists of smooth functions (3, going 
rapidly to zero when |x| — )• oo such that (3{x) G ^si^) &i^d P{(3) G Ts{Q). The last means that 

i7Tsd^(3 - iS.,G-^d-\d^l5{x), S(x)) G f^(x), 

that is, for arbitrary smooth function X(x) such that X G and going rapidly to zero when 
|x| — )• oo we have 

(ivr^a^/^ - iS^G-^d-\d^p{x),X{x)) = 0. (34) 

The space T-^{Q)s consists of smooth (3{x) such that /3 G f)^^, going rapidly to zero when — )• oo 
and satisfying 

{{iTTsd.P -iS^G-'d-\d.,P{x),S{x)),Y{x))) =0 

for each smooth function Y{x), Y G fg{x) going rapidly to zero when |x| — )• oo. Arguments 
similar those used to prove Haar's lemma in the variational calculus show that for each x 

{i7rsd,(3 - iS,G-'d-'{d,(3{x),S{x)),Y{x)) = 0, (35) 

where Y{x) is as above. Then, if /3 G X*{P)s H T^{Q)s, we shall have simultaneously (34) 
and (35) so /3 G kerP^. Thus the first condition of the first restriction theorem is fulfilled. In 
order to see that the second condition also holds, we introduce 
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Lemma 1. The operator P has the properties 

Psmfs))cd{fs), Psmfs))cm's)- 

The proof of the lemma is obtained easily since the spaces and are invariant with respect 
to vr^. 

Using the lemma, suppose /3 G f);^ is a smooth function going rapidly to zero when |x| — )■ oo. 
Clearly, we can write it uniquely into the form 

/3 = /3o + /3i, (3oems), f^i^ms)- 

Then PsiPo) G 5(f^), Ps{(3i) G d{fs) and we see that Ps/^o G X*{P)s, PsPi G T^{Q)s- So the 
second requirement of the first restriction theorem is also satisfied and P allows restriction. If 7 
is a 1-form on Q, that is 7 G 5(f5)i the restriction P' is given by 

P'{^) = msd^^-\Si,^{^G)-^d-\d^^{x),Si{x)). (36) 

We summarize these facts into the following 

Theorem 5. The manifold of potentials Q is endowed with a P-N structure, defined by the 
Poisson tensor P' (36) and the restriction of the Nijenhuis tensor N'^ . Explicitly, the restriction 
of iV^ is given by NqNi where: 

No{X) = iTTsd^iads^ X) - iSi.^{'Gr^d~^{d.,{ads' X), Si(x)), 

Xedifs), No{X)ed{fs), 
Ni{X) = insd^iads^ X) - iS,d~\d,{ads^ X),Soix)), 

Xedifs), No{X)ems)- 

Since (Nq)* = A^, (A'^i)* = A5j_ the above theorem gives a geometric interpretation of the 
operators A^'^. 

Finally, let us calculate the Poisson bracket for the Poisson structure P' . If Hi, H2 are two 
allowed functionals and we introduce the row (Si;^;, ^^^ix) with elements {Sr-x, ^^)ix)i r < n, 
r = 1, 3, 5, . . . , we shall obtain 



{Hi,H2} = i((dx 



5Hi\ 6H2 



5S r 6S 



+ 2 



dx. 



4.3 The hierarchies of integrable equations 

The geometric theory is incomplete without giving (at least partially) the hierarchies of equations 
(vector fields) that correspond to the P-N structure described in the above theorem, that is the 
fields that are fundamental for the P-N structure. In order to discuss this issue let us start with 
the fundamental fields for the P-N structure in the general case (without reductions). As it is 
known, see [10], the fields of the type Xh{S) = [H,S], H G I) are fundamental for the P-N 
structure. This is because the original tensors P, Q (see (19)) are covariant with respect to the 
one-parametric group of transformations 

(5) = Ad(expt//)5, HGi), tGM, 

and the submanifold Af is invariant with respect to tpj^ . The one-parametric group (p^ corre- 
sponds to the vector fields Xh which are tangent to M. According to the theory of the P-N 
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manifolds, if X is a fundamental field then all the fields N^X, p = 1,2, . . . are also fundamental 
and commute, that is they have zero Lie brackets. Thus to each € f) corresponds a hierarchy of 
fundamental fields N^[H, S] and we obtain n — 1 independent families of fundamental fields. For 
different p and different G f) the fields from these families also commute. As a matter of fact 
this is commonly referred to as 'the geometry' behind the properties of the hierarchies (8). We 
can cast these hierarchies in a different form. Using the expression (21) and taking into account 
that ad^^ ad5 H = 7rs{H) after some simple transformations in which we use the properties of 
the Killing form and Proposition 1 we get 

N{XH)iS) = iTTsd^iirsiH)) - [S,G~^d^\d,i7rsiX)),Six)) 

= -iS^.G-^(F,S(±oo)) = -S^G-^(-H',iS(±oo)). (37) 

The last expression shows that N{Xh){S) is a linear combination of the fields Si-x = Sx, S2-x = 
{S'^)x, . . • , Sn-i;x = {S^~^)x- Since the evolution equations corresponding to the fields Xh are 
linear one can say that the hierarchies of the fundamental fields that correspond to nonlinear 
equations are generated by the above vector fields. Now, let us assume that one has reductions 
as in the above and let us first consider what happens when we restrict to isu(n). In that case 
we must take H S su(n), that is, iH must be a diagonal, traceless real matrix. Now we can see 
that the coefficients in front of the fields S^-x are real since su(n) is a compact real form of sl(n) 
and on it the Killing form is real so both G and (i?, iS(iboo)) are real. 

When we have the additional restriction defined by the automorphism h the situation is 
more complicated. Since h(S) = —S and h{H) = H for any H £ I) we have h^Xn) = —Xh so 
Xh e ^{^^s)- Thus N{Xh) G ^{%s) and more generally for p = 0, 1, 2, . . . 

N^^Xh = {NoN.yXH e i?(i(i^), N^p+'Xh = {NiNoTNiXh G 5(%^). (38) 

Remark 11. The fact that N{Xh) G di°^s) 

can also be seen easily from (37). Indeed, because 
H is an element of the Cartan subalgebra it is orthogonal to all elements of fi and in (37) we 
have a linear combination of the elements Sm with odd index m, that is —SxG''^{H, iS(iboo)) = 
-Si;,iG-i(i7,iSi(±oo)). 

Thus when we have a reduction defined by the automorphism h the first of the series in (38) 
consists of vector fields tangent to the manifold Q while the vector fields from the second series 
are not tangent to Q. 

Naturally, the vector fields from the first series are the candidates for being fundamental 
fields of the P-N structure on Q and, discarding the first fields in the hierarchies, we see that 
we have fields of the type 

N^''S2j-i-x = {NoNirS2j-i-x- 

These are not, however, all the fundamental fields that we can produce. The fields of the 
type NS2j-x are fundamental for N and hence are fundamental for N'^ also. Besides, since 
S2j;x S 5^('^fl5) we have NS2j-x £ i?(^^)5)- As easily checked, the fields NS2j-x are tangent to Q. 

Thus, finally, the hierarchies of the the fundamental fields of the P-N structure when we 
restrict to Q are generated by the fields 

N^'''S2j-i;x = {NoNi)PS2j-i;x, N^''+^S2i-x = {NoNi)PNoS2i;x, 
l<j<k, l<l<s, p = 0,1,2,..., 

where k and s are such that 2A; — 1 is the largest odd number less than n and 2s is the largest 
even number less than n. These numbers define the sizes of the Gram matrices ^G and ^G. For 
example, in the case of the algebra sl(3,C), we have the fields 

N^PSx = {NoN,rSx, N^P+Hs\ = {NoNirNoiS\, p = 0,1,2.... (39) 
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These fundamental fields give rise to hierarchies of integrable equations of the type St = F{S) 
where F{S) is a finite linear combination of the fundamental fields (39), see [7]. 

5 Conclusions 

We have been able to show that the geometric interpretation known for the the recursion opera- 
tors related to the generalized Zakharov-Shabat system on the algebra sl(n, C) in pole gauge 
holds also in the presence of Z2 x Z2 reductions of certain classes and we have explicitly calculated 
the recursion operators. It is an interesting question whether analogous results could be obtained 
for Zakharov-Shabat type systems on the other classical Lie algebras. 
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